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Abstract 

In this article, we study the structures of the pseudoscalar mesons n, K 
and the scalar diquarks U a , D a , S a in the framework of the coupled rainbow 
Schwinger-Dyson equation and ladder Bethe-Salpeter equation with the con- 
fining effective potential. The u, d, s quarks have small current masses, and 
the renormalization is very large, the mass poles in the timelike region are 
absent which implements confinement naturally. The Bethe-Salpeter wave- 
functions of the pseudoscalar mesons tt, K and the scalar diquarks U a , D a , S a 
have the same type (Gaussian type) momentum dependence, center around 
zero momentum and extend to the energy scale about q 2 = lGeV 2 which 
happen to be the energy scale for the chiral symmetry breaking, the strong 
interactions in the infrared region result in bound (or quasi-bound) states. 
The numerical results for the masses and decay constants of the n, K mesons 
can reproduce the experimental values, the ground state masses of the scalar 
diquarks U a , D a , S a are consistent with the existing theoretical calculations. 
We suggest a new Lagrangian which may explain the uncertainty of the masses 
of the scalar diquarks . 
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1 Introduction 

The discovery of the pentaquark state + (154O) has opened a new field of strong 
interaction and provides a new opportunity for a deeper understanding of the low 
energy QCD. Intense theoretical studies have been motivated to clarify the quantum 
numbers and to understand the under-structures of the pentaquark state G + (1540) 
p. Although the existence of the + (154O) is uncertain and is a subject of contro- 
versy now, the 6 + (1540) has already contributed to hadron spectroscopy. Unlike 
the chiral soliton model, the quark models take the constituent quarks or quark 
clusters as the elementary degrees of freedom, there exist a great number of possible 
quark configurations satisfy the Fermi statistics and the color singlet condition for 
the substructures of the pentaquark state + (154O) if we release stringent dynamical 
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constraints. In fact, the multiquark states are many-body problems, they are very 
difficult to solve. Whether or not the quarks can cluster together to form diquarks 
is of great importance theoretically, if we take the diquarks as the basic constituents 
(here the "basic constituents" does not mean they are asymptotic states, they just 
exist inside the baryons or multiquark states with typical length), the problems will 
be greatly simplified, for example, the \ octet and | decuplet baryons can be 
taken as quark-diquark bound states in the Faddeev approximation [2|, the nonet 
scalar mesons below lGeV can be taken as 4-quark states (qq)s(qq)3 with scalar di- 
quarks or pseudoscalar diquarks as their basic constituents [3]; furthermore, we can 
obtain more insight into the relevant degrees of freedom and deepen our understand- 
ing about the underlying dynamics that determines the properties of the baryons 
and exotic multiquark states. A typical quark model for the pentaquark states is the 
Jaffe-Wilczek's diquark-diquark-antiquark model In this model, the scalar di- 
quarks U a = e abc dl(x)C-f 5 s c (x), D a = e abc ul(x)C-f 5 s c (x) , S a = e abc ul(x)C^d c (x) 
are taken as the basic constituents. They belong to the antitriplet 3 representation 
of both the color SU(3) C group and flavor SU(3)j group, in the color superconduc- 
tivity theory, the attractive interactions in this channel lead to the formulation of 
nonzero condensates and breaking of both the color and flavor 577(3) symmetries 
j^]. The scalar diquarks correspond to the 1 So states of the diquark systems, the 
one-gluon exchange force and the instanton induced force can lead to significant 
attractions between the quarks in the + channels jH]. The pseudoscalar diquarks 
do not have nonrelativistic limit, can be taken as the 3 Pq states. As the instanton 
induced force results in strong attractions in the scalar diquark channel and strong 
repulsions in the pseudoscalar diquark channel, if the effects of the instanton are 
manifested, we prefer the S type diquark to the P type diquark in constructing 
interpolating currents in the QCD sum rules 

In this article, we take the point of view that the scalar diquarks are quasi- 
bound states of quark-quark system and study the ground state mass spectrum 
within the framework of the coupled Schwinger- Dyson equation (SDE) and Bethe- 
Salpeter equation (BSE). The coupled rainbow SDE and ladder BSE have given a lot 
of successful descriptions of the long distance properties of the low energy QCD and 
the QCD vacuum (for reviews, one can see Refs.[7|l8ll§llT0]). The SDE can naturally 
embody the dynamical symmetry breaking and confinement which are two crucial 
features of QCD, although they correspond to two very different energy scales jit; 
IT2] . On the other hand, the BSE is a conventional approach in dealing with the two- 
body relativistic bound state problems |T3] • From the solutions of the BSE, we can 
obtain useful information about the under-structures of the mesons and diquarks, 
and obtain powerful tests for the quark theory. However, the obviously drawback 
may be the model dependent kernels for the gluon two-point Green's function and 
the truncations for the coupled divergent SDE and BSE series in one or the other 
ways [Tl]. Many analytical and numerical calculations indicate that the coupled 
rainbow SDE and ladder BSE with phenomenological potential models can give 
model independent results and satisfactory values CHI - The usually used effective 
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potential models are confining Dirac 5 function potential, Gaussian distribution 
potential and flat bottom potential (FBP) [T71ll8| IT^]. The FBP is a sum of Yukawa 
potentials, which not only satisfies chiral invariance and fully relativistic covariance, 
but also suppresses the singular point that the Yukawa potential has. It works well 
in understanding the dynamical chiral symmetry breaking, confinement and the 
QCD vacuum as well as the meson structures, such as electromagnetic form factors, 
radius, decay constants [HI I2HI EH 122]- In this article, we use the FBP to study 
the ground state mass spectrum of the scalar diquarks without fine tuning such as 
modifying the infrared behavior for the heavy quark systems. 

The article is arranged as follows: we introduce the FBP in section II; in section 
III, IV and V, we solve the rainbow SDE and ladder BSE, explore the analyticity of 
the quark propagators, study the dynamical symmetry breaking and confinement, 
finally obtain the mass spectrum of the 7r, K mesons and the scalar U a , D a , S a 
diquarks, and the decay constants of the tt, K mesons; section VI is reserved for 
conclusion. 



2 Flat Bottom Potential 

The present techniques in QCD calculation can not give satisfactory large r behavior 
for the gluon two-point Green's function to implement the linear potential confine- 
ment mechanism, in practical calculation, the phenomenological effective potential 
models always do the work. The FBP is a sum of Yukawa potentials which is an 
analogy to the exchange of a series of particles and ghosts with different masses 
(Euclidean Form), 

n 

where N stands for the minimum value of the masses, p is their mass difference, and 
a,j is their relative coupling constant. Due to the particular condition we take for 
the FBP, there is no divergence in solving the SDE. In its three dimensional form, 
the FBP takes the following form: 

JL e -(N+jp)r 

V(r) = -$>^— . (2) 

j=0 

In order to suppress the singular point at r = 0, we take the following conditions: 



V(0) = constant 
dV(0) d 2 V(0) d n V(0) 



dr dr 2 dr n 



0. (3) 



The cij can be determined by solving the equations inferred from the flat bottom 
condition in Eq.(3). As in previous literature [IU EH 123 HU 122] , u is set to be 9. 
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3 Schwinger-Dyson equation 



The SDE can provide a natural framework for studying the nonperturbative prop- 
erties of the quark and gluon Green's functions. By studying the evolution behavior 
and analytic structure of the dressed quark propagators, we can obtain valuable in- 
formation about the dynamical chiral symmetry breaking and confinement. In the 
following, we write down the rainbow SDE for the quark propagator, 

/d 4 k X a X a 
-^^l^ s ( k )lu^G^{k-p), (4) 

where 

S-\p) = iA(p 2 )~f -p + B(p 2 ) = A(p 2 ){i~f -p + m(p 2 )], (5) 
G, v {k) = (5, u -^)G(k 2 ), (6) 

and rh u ^ yS stands for the current quark mass that breaks chiral symmetry explicitly. 
For a short discussion about the full SDE for the quark propagator, one can consult 
Ref.jUl. 

In this article, we assume that a Wick rotation to Euclidean variables is allowed, 
and perform a rotation analytically continuing p and k into the Euclidean region. 
Alternatively, one can derive the SDE from the Euclidean path-integral formulation 
of the theory, thus avoiding possible difficulties in performing the Wick rotation 
[23] . As far as only numerical results are concerned, the two procedures are equal. 
In fact, the analytical structures of quark propagators have interesting information 
about confinement, we will make detailed discussion about the propagators of the 
u, d and s quarks in section V. 



4 Bethe-Salpeter equation 

The BSE is a conventional approach in dealing with the two-body relativistic bound 
state problems [T3]. The precise knowledge about the quark structures of the mesons 
and diquarks can result in better understanding of their properties. In the following, 
we write down the ladder BSE for the scalar diquark quasi-bound states with two 
quarks of different flavor j2H , 

r s (g,P)= (7) 

Here T denotes matrix transpose, the S(k) is the quark propagator, G^ik) is the 
gluon propagator, P M is the four-momentum of the center of mass of the scalar 
diquark, is the relative four-momentum between the two quarks, 7^ is the bare 
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quark-gluon vertex, and T^(q, P) is the Bethe-Salpeter amplitude of the quasi-bound 
state (or diquark). The £ is the center of mass parameter which can be chosen to 
vary between and 1, for the S a diquark, £ = |, for the U a and D a diquarks, as 
the current quark masses m s > m u and m s > ma, £ is about \. We can introduce 
an auxiliary amplitude T^(q,P) to facilitate the calculation, 

iffo P) = P) C, (8) 

here C = 7274 is the charge conjugation matrix. The auxiliary amplitude Ig (g, P) 
satisfies the following equation, 

If (?, P) = ~ y J J^y G ^ - *)7, S(k + £P) if (*, P)S(* - (1 - Z)P) lv . (9) 

It is obviously the above equation is identical to the BSE for the pseudoscalar 
mesons but a reduction in the coupling strength, 3 — > |. We can introduce the 
Bethe-Salpeter wavefunction (BSW) Xgg for the quasi-bound states, 

Xqq (q, P) = S(q + £P)lf (p, P)5(? - (1 - OP), (10) 
to relate with our previously works on the pseudoscalar meson's BSE, 

87T f H 4 k 

S-\q + ZP)x qq (q,P)S-\q-(l-OP) = -Y J (2^X W (*, P )l^, v {q - k), (11) 

1fi7r /" r/ 4 A- 

S-\q + £P)x 9g -(g, P)^? - (1 - OP) = ~— J J^y^X qq (k, P)l v G» v {q - fc),(12) 

We can perform the Wick rotation analytically and continue q and /c into the 
Euclidean region 2 . The BSWs of the scalar diquarks x qq an d pseudoscalar mesons 
X qq have the same Dirac structures, can be signed by the notation x(g, P). In the 
lowest order approximation, the BSW x(<2S P) can be written as 

X (q, P) = 75 [iF?(q, P) + 7 • PP 2 °(<7, P) + 7 • qq • PFafo P) + *[7 • 5, 7 • P]P 4 °(?, ^)] -(13) 

In solving the BSEs, it is important to translate the wavefunctions P™ into the same 
dimension, 

P° -> A°P 1 °, P 2 ° - A 1 ^ , P3 1 - A 3 P 3 \ P 4 ° - A 2 P 4 °, q -> g/A, P^P/A, 

here the A is some quantity of the dimension of mass. The ladder BSEs for the 
scalar diquarks and pseudoscalar mesons can be projected into the following four 
coupled integral equations, 

Y,H(i,j)Ff>\q,P) = W^X(m), (14) 

j 3 J 

2 To avoid possible difficulties in performing the Wick rotation, one can derive the BSE from 
the Euclidean path-integral formulation of the theory. 
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the expressions of the H(i,j) and K(i,j) are cumbersome and neglected here. 

We can introduce a parameter A(P 2 ) and solve the above equations as an eigen- 
value problem. If there really exist a quasi-bound state of two-quark, the masses of 
the diquarks can be determined by the condition A(P 2 = — M 2 ) = 1, 

TT,H{i,j)lf\q,P) = \(P 2 )J2 fd*kK(i,j). (15) 

3 3 J 

Here we will take a short digression and give some explanations for the ex- 
pressions of H(i,j) . The H(i,jys are functions of the quark's Schwinger-Dyson 
functions (SDF) A(q 2 + £ 2 P 2 + 2£q ■ P) ,B(q 2 + £ 2 P 2 + 2£q ■ P) , A(q 2 + (1 - £) 2 P 2 - 
2(1 - £)q ■ P) and B(q 2 + (1 - £) 2 P 2 - 2(1 - g)q ■ P) . The relative four- momentum 
q is a quantity in the Euclidean spacetime while the center of mass four-momentum 
P must be continued to the Minkowski spacetime i.e. P 2 = — m 2 K ua Da Sa , this 
results in that the q ■ P varies throughout a complex domain. It is inconvenient to 
solve the SDE with the resulting complex values of the quark momentum. We can 
expand the A and B in terms of Taylor series of q ■ P, for example, 

A(q 2 + £ 2 P 2 + £q ■ P) = A{q 2 + £ 2 P 2 ) + 2£A(q 2 + ?P 2 )'q ■ P + ■ ■ ■ . 

The other problem is that we can not solve the SDE in the timelike region as the 
two-point gluon Green's function can not be exactly inferred from the SU (3) color 
gauge theory even in the low energy spacelike region. In practical calculations, we 
can extrapolate the values of the A and B from the spacelike region smoothly to the 
timelike region with suitable polynomial functions. To avoid possible violation with 
confinement in sense of the appearance of pole masses q 2 = —m 2 (q 2 ) in the timelike 
region, we must be care in choosing the polynomial functions 18 . 

Finally we write down the normalization condition for the BSWs of the pseu- 
doscalar mesons, 

Nc J j^ Tr v^ x ^ p)szl + ^ lx ^ p) ^s = 2P ^ (16) 

here \ — 74X + T4j S+ — S(q + CP) an d S- = S(q — (1 — £)P). In this article, the 
parameters of FBP are fitted to give the correct masses and decay constants for the 
pseudoscalar mesons, 7r and K, the normalization condition is needed. 

5 Coupled rainbow SDE and ladder BSE, and the 
mass spectrum 

In this section, we study the coupled equations of the rainbow SDE and ladder 
BSE for the pseudoscalar mesons (tt and K) and scalar diquarks (U a , D a and S a ) 
numerically, the final results for the SDFs and BSWs can be plotted as functions of 
the square momentum q 2 . 
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In order to demonstrate the confinement of quarks, we have to study the analyt- 
icity of the SDFs of the u, d and s quarks, and prove that there are no mass poles 
on the real timelike q 2 axial. In the following, we take the Fourier transform with 
respect to the Euclidean time T for the scalar part (S s ) of the quark propagator 



S s( T ) - / TTT etq4T oAO,o\ ; 5?i73Fs l"?=o> ( 17 ) 



2tt g 2 A 2 (g 2 ) + 5 2 (g 2 )' T=0 ' 

where the 3- vector part of q is set to zero. If S(q) has a mass pole at q 2 = —m 2 (q 2 ) 
in the real timelike region, the Fourier transformed S*(T) would fall off as e~ mT for 
large T or log S* = —mT. In our numerical calculations, for small T, the values 
of S* are positive and decrease rapidly to zero and beyond with the increase of T, 
which are compatible with the result (curve tendency with respect to T) from lattice 
simulations j2Zj ; for large T, the values of S* are negative, except occasionally a 
very small fraction of positive values. The negative values for S* indicate an explicit 
violation of the axiom of reflection positivity [2B|> m other words, the quarks are 
not physical observable i.e. confinement. As colored quantity, the diquarks should 
also be confined and not appear as asymptotic states; the demonstration of their 
confinement is beyond the present work. 

The u, d and s quarks have small current masses, the dressing or renormaliza- 
tion is very large and the curves of the SDFs are steep, which are corresponding 
to the dynamical chiral symmetry breaking phenomenon for the light quarks. At 
zero momentum, m u (0) = m^(0) = 0A54GeV and m s (0) = 0.68AGeV , which are 
compatible with the constituent quark masses in the literature. From the solutions 
of BSEs for the tt, K mesons and U a , D a , S a diquarks as eigenvalue problems, we 
can obtain the masses for those pseudoscalar mesons and scalar diquarks, 

M n = 135MeV, M K = A98MeV, 
M S a = 0.76GeV, M ua = 0.98GeV, M D a = 0.98MeV. (18) 

It is obviously 

M s « < m u {0) + m d (0), M Ua < m d {0) + m s (0), M D * < m u {0) + m,(0), (19) 

M ua - M s * w m s {0) - m u (0) w 0.22GeV. (20) 

The attractive interaction between the quarks in the color and flavor SU(3) 3 channel 
can lead to the quasi-bound states in the infrared region. The appearance of the 
diquarks is closely related to the dynamical symmetry breaking phenomenon, the 
mass splitting among the U a , D a and S a diquarks originate from the mass splitting 
among u, d and s quarks. 

The existing theoretical calculations for the masses of the scalar diquarks vary in 
a large range, M qq = (0.4 - 0.7)GeV with QCD sum rules |2H|; M qq ~ 0.5GeV with 
random instanton liquid model jSHj; M qq = (0.42±0.03)GeV with random instanton 
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Figure 1: Un- normalized BSWs for the it meson. 



liquid model (HH, M qq = 0.234GW with Nambu-Jona-Lasinio Model [221; M s a = 
0.7AGeV,Mua = M D a = 0.88GeV with BSE j2£J; M 5 a = 0.82GeV, Mj/a = M D a = 
l.lOGeV with BSE ; M qq = 0.692GeV with global color model [M|; 0.14 < 
M qq < 0.74GeV with assumption of mass functions M gg « 0.7GeV with lattice 
QCD [35]. There are large uncertainties for the masses of the quasi-bound states, 
U a , D a , and S a . As colored quantities, the diquarks may have gauge interactions 
with the gluon field as fundamental scalar field with the following Lagrangian, 



field. The direct color interactions between the scalar diquarks and the gluons may 
modify the mass M qq significantly. 

From the plotted BSWs (see Fig.l for the 7r meson and Fig. 2 for the 5* a diquark as 
examples), we can see that the BSWs of the pseudoscalar mesons and scalar diquarks 
have the same type (Gaussian type) momentum dependence while the quantitative 
values are different from each other. Just like the qq , qQ and QQ pseudoscalar 
mesons [HJ |20], the gaussian type BSWs of the scalar diquarks center around zero 
momentum and extend to the energy scale about q 2 = lGeV 2 which happen to be 
the energy scale for the chiral symmetry breaking, the strong interactions in the 
infrared region result in quasi-bound states, U a , D a and S a . The BSWs of the n 
and K mesons can give satisfactory values for the decay constants which are defined 
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Figure 2: Un-normalized BSWs for the S a diquark . 



by 

= N c Jtt[ WX %P)]^ (23) 

here we use 7r to represent the pseudoscalar mesons, 

U = !27MeV; f K = 161MeV. (24) 

There are negative voice for the existence of the scalar diquark states in the 
infrared region [36j. The coupled rainbow SDE and ladder BSE are particularly 
suitable for studying the flavor octet pseudoscalar mesons and vector mesons, the 
next-to-leading-order (NLO) contributions from the quark-gluon vertex have a sig- 
nificant amount of cancellation between repulsive and attractive corrections. How- 
ever, for the diquarks and scalar mesons, the large repulsive corrections from the 
NLO contributions can significantly change the scalar meson masses and correspond- 
ing BSWs; the diquark quasi-bound states found in the ladder BSE will disappear 
from the spectrum. In this article, we do not take the point of view that the scalar 
diquarks exist in the strong interaction spectrum as asymptotic states, and our con- 
clusion do not conflict with Ref . |36j : the confinement precludes the observation of 
the free colored diquarks, the quark-quark can correlate with each other in the color 
and flavor 3 channels inside the baryons and multi-quark states with typical length 
I = -tj— . In fact, the NLO contributions from the quark-gluon vertex are extremely 
difficult to take into account if we go beyond the infrared dominated S function ap- 
proximation for the gluon kernel, which is obviously violate the lorentz invariance; 
the lorentz invariant and model-independent treatments still lack in the literatures. 
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In calculation, the values of current quark masses are taken as rh u = rh d = 
QMeV and rh s = 150MeV; the input parameters for the FBP are N = 1.0A, 
V(0) = -17.0A, p = 6.0A and A = 200MeV, which are fitted to give the correct 
masses of the 7r and K mesons. In this article, we deal with only the light flavor 
quarks, the FBP can give satisfactory results without fine tuning, such as modifying 
the infrared behavior for the heavy quarks c and b. 

6 Conclusion 

In this article, we study the under-structures of the pseudoscalar mesons n , K 
and scalar diquarks U a , D a , S a in the framework of the coupled rainbow SDE and 
ladder BSE with the confining effective potential (FBP). After we solve the coupled 
rainbow SDE and ladder BSE numerically, we obtain the SDFs and BSWs of the 
pseudoscalar mesons n , K and scalar diquarks U a , D a , S a , and the corresponding 
ground state mass spectrum. The u, d and s quarks have small current masses, the 
dressing or renormalization for the SDFs is very large and the curves are steep which 
indicate the dynamical chiral symmetry breaking phenomenon for the light quarks 
explicitly. The mass poles in the timelike region are absent which implement the 
confinement naturally. The BSWs of the pseudoscalar mesons and scalar diquarks 
have the same type (Gaussian type) momentum dependence while the quantitative 
values are different from each other. The gaussian type BSWs center around zero 
momentum and extend to the energy scale about q 2 = lGeV 2 which happen to be the 
energy scale for the chiral symmetry breaking, the strong interactions in the infrared 
region result in bound (or quasi-bound) states. Our numerical results for the masses 
and decay constants of the n, K mesons can reproduce the experimental values, 
the ground state mass spectrum of scalar diquarks are consistent with the existing 
theoretical calculations. The mass splitting among the U a , D a and S a diquarks 
originate from the mass splitting among u, d and s quarks. M$a < m u (0) + m^(0), 
M\ja < m^O) + m s (0) , Mi)a < m u (0) + m s (0) , the attractive interaction between 
the color and flavor SU (3) 3 quarks can lead to the quasi-bound states in the infrared 
region. Once the satisfactory SDFs and BSWs of the scalar diquarks are known, we 
can use them to study a lot of important quantities involving the multiquark states. 
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